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Abstract. To investigate the stability of protoneutron stars in their early evolution, the minimum gravitational mass plays a 
(2JT)| fundamental role. This quantity depends upon the temperature profile assumed. We study within a static approach the stability of 

a protoneutron star. In particular we focus on a suitable temperature profile suggested by dynamical calculations. We consider a 
protoneutron star as composed of an isothermal core and an isentropic outer part. To describe physical properties of the interior 
we employ a microscopically derived equation of state for nuclear matter. For the outer part we employ the Lattimer-Swesty 
equation of state. The global structure is studied. The assumed temperature profile turns out to give a range of stability which 
supports temperature values in line with those coming from dynamical calculations. The maximum mass instead depends only 
upon the equation of state employed. 
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O 1- Introduction 

\Q . 

• Soon after a supernova explosion a protoneutron star (PNS) is formed and constitutes for several tens of seconds a transitional 
| state to either a neutron star (NS) or a black hole (BH). In the early stage of evolution a PNS is dominated by neutrino diffusion 

-i— > . | |Burrows & Lattimer (198 6) I from the core to the surface, resulting firstly in deleptonization, which heats the star generating 
i— h ' temperature values ranging from 30 to 50 MeV. 

y . In this note we will focus on the study of the structure of a PNS in its early stage in the framework of a static approach. We 
assume, prompted by several dynamical calculations ( |Burrows & Lattimer (1986)1 |Pons et al.(1999)| |Sumiyoshi et al.(2005)^ , 

• • . that in this stage the neutrinos are trapped and the temperature is high. In particular we will study those configurations where the 
. ^ ' core has a temperature greater than 10 MeV and the outer part of the star is characterized by a high value of entropy per baryon. 

• When one deals with a static description of these configurations a problem occurs: what temperature profile shall we assume? 
5— i ' This key question will be the central point of this note. 

. 5^ i The basic input for each description of PNS is the nuclear matter equation of state (EoS). For the neutrino-rich core the 
non-relativistic Brueckner-Hartree-Fock (BHF) EoS at finite temperature (Baldo(1999)jl was implemented. We will not give 
here further details on the EoS implemented and we refer the reader to the original works (Baldo and Ferreira(1999)l. For the 
high-entropy envelope we implement the EoS developed by Lattimer-Swest y (LS) (|Latt imer & S westy(T997)| . A description of 
PNS within the BHF approach has already been given in a preceding work tNico traeTaL(2006)t . In this sense this note can be 
considered as a refinement of that work with focus on a particular choice for the temperature and entropy profiles. 
The note is organized as follows. In section[2]we present a new static approach for PNS which is built taking care as much as 
possible of all information coming from dynamical calculations mentioned above. Sectionals devoted to discussing the results 
implied by the approach. Finally in section@]we draw our conclusions. 



2. Model 

Calculations of static models of protoneutron stars should be considered as a first step to describe these objects. In principle 
the temperature profile has to be determined via dynamical calculations taking into account neutrino transport properly. Many 
static approaches have been developed in the past decade ( Praka shet al.(1997)| , implementing several finite temperature EoS and 
assuming an isentropic or an isothermal ( Pra kash et al.(199^)| l profile throughout the star. Nevertheless these simple assumptions 
on the temperature profile lead to some problems. For the case of an isentropic profile static models can predict, in pure nucle- 
onic PNS, values of central temperature higher than those predicted by dynamical simulations, or they are unable to reproduce 



2 O. E. Nicotra: Structure of protoneutron stars within a static approach 



40 



30 



20 



10 





1 






- 


v i 


1 1 1 1 " 


/- 








0.8 

0.6 .' "7 


";= 
















- ^ 


















1 1 






0,002 


3.004 0.006 0,008 0.01 -' / / 






p[fm] 






s=2 








s=4 








— - s=6 








s=8 








s=10 
















1 


1 





I — — — — 

le-06 le-05 0,0001 0,001 0,01 



p[fm ] 



Fig. 1. Temperature profile for different fixed values of entropy per baryon s as a function of the baryon density for the Lattimer- 
Swesty EoS. The corresponding pressure as a function of baryon density is reported in the inner panel. 



correctly minimum masses for the isothermal case. For the latter case much better results were obtained using a more realistic 
prescription (Gondek et al.(1997)i, which also allows to predict location and temperature of the neutrinosphere. A very interest- 
ing assumption on temperature and entropy profiles can be found also in |Strobel et al.(1999)| where the entire PNS evolution 
consists of four stages. 

In our model we assume that a PNS is composed of a hot, neutrino opaque, and isothermal core separated from an outer cold 
crust by an isentropic, neutrino-free intermediate layer, which will be called the envelope throughout the paper. 



2.1. isothermal core 

For a PNS core in which the strongly interacting particles are only baryons, its composition is determined by requirements of 
charge neutrality and equilibrium under weak semileptonic processes, 

B x -^B 2 + l + v,, B 2 + /->fii +v, , (1) 

where B\ and Bi are baryons and / is a lepton (either an electron or a muon). Under the condition that neutrinos are trapped in 
the system, the beta equilibrium equations read explicitly 

yu„ - fi p = fx e - fX Ye =/J. M - /iv„ ■ (2) 

Because of trapping, the numbers of leptons per baryon of each flavour {I = e,/j.), Yi = xi - Xj + x v , - x v are conserved. 
Gravitational collapse calculations of the iron core of massive stars indicate that, at the onset of trapping, the electron lepton 
number is Y e ^ 0.4; since no muons are present at this stage we can impose also = 0. For neutrino free matter we just set 
fj. Vl — in eq.@ and neglect the above constraints on lepton numbers. We assume a constant value of temperature throughout 
the core and perform some calculations for a value of temperature ranging from to 50 MeV, with and without neutrinos. The 
EoS employed is that of the BHF approach at finite temperature. Plots for the chemical composition and pressure at increasing 
density and temperature can be found in [Nicotra et al.( 2006) 

2.2. isentropic envelope 

The condition of isothermality adopted for the core cannot be extended to the outer part of the star. Dynamical calculations suggest 
that the temperature drops rapidly to zero at the surface of the star; this is due to the fast cooling of the outer part of the PNS where 
the stellar matter is transparent to neutrinos. Moreover, in the early stage, the outer part of a PNS is characterized by a high value of 
the entropy per baryon. Essentially, during the first 10 seconds, the entropy profile decreases from the surface to the core starting 
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Fig. 2. The gravitational mass Mq versus central energy density e c (left panel) and versus stellar radius R (right panel) for several 
PNS static configurations (see legend) as discussed in the paper. The gravitational mass is in units of solar masses Mq, central 
energy density in units of eo (energy density at saturation), and R in Km. 

from values of 6 H- 10 in units of Boltzmann's constant ( Burrow^^I^attim er (1986)||Pons et al.(1999)||Sumiyoshi et al.(2005)) . 
For a low enough core temperature (T < 10 MeV) in |Nicotn ^H^2 006) | a tem perature profile in the shape of a step function 
was assumed, joining the hadronic EoS (BHF) with the BPS i Bayn^i al.(1971)) plus FMT ( |Feynman et al.(19~9 )it EoS for the 
cold crust. This approach does not work anymore if the core temperature T core is greater than 10 MeV, leading to a value for the 
minimum gravitational mass which rises too fast with increasing T core . This actually makes the PNS already unstable at T core 
slightly above 30 MeV in the neutrino-free case. 

In order to avoid this problem we consider an isentropic envelope in the range of baryon density from 10~ 6 frrr 3 to 0.01 fm" 3 
based on the EoS of LS ( |Lattimer & Swesty(1997)| l with the incompressibility modulus of symmetric nuclear matter K - 220 
MeV (LS220). Fixing the entropy per baryon s to the above values we get the temperature profiles reported in FigQ they 
explicitly suggest a natural correspondance between the entropy of the envelope s env and T core . Moreover, for the curves displayed 
in Fig^beta equilibrium is imposed assuming a neutrino-free regime. We thus have the EoS of the envelope with a temperature 
profile which rises quickly to typical values of T core , as the density increases, in a narrow range of baryon density, just fixing the 
entropy per baryon at a value consistent with dynamical calculations. 

Now we need to join together this EoS and that of the star interior. Since the energy of neutrinos, emerging from the interior, 
possesses some spreading ( Bethe ( 1 990) I, we do not think that the location of the neutrino-sphere, affected by the same spreading, 
is a good criterion to fix the matching density between core and envelope. Neutrino transport properties indeed vary quite a lot 
during the PNS evolution, they depend on the stellar matter composition and change with the leptonic flavour leading to different 
locations of the energy spheres and transport spheres for different neutrino types (Janka (1995) I. In addition, the temperature 
drop regards also those profiles chosen as initial conditions for dynamical calculations and borrowed from results on iron core 
collapse simulations, thus establishing its independence from the subsequent neutrino diffusion. All this suggests the possibility 
to have a blurring region inside the star where we are free to choose the matching density. Conversely, if we consider s em , as a 
free parameter in order to join core and envelope energy and free energy densities, we recognize soon that only small changes 
around p em , ^ 0.01 firT 3 are allowed in order to maintain s env consistent with the dynamical expectations. Nevertheless, to build 
up our model of PNS, a fine tuning of s em , around the values suggested by Figgis performed in order to have an exact matching 
between core and envelope of all the thermodynamic quantities (energy density, free energy density, and temperature). In this 
sense we consider s env as a free parameter. 

3. Results 



The stable configurations of a PNS can be obtained from the hydrostatic equilibrium equation of Tolman, Oppenheimer, and 
Volkov (Shapir (7&Teukolsky(l 983)| |Zeldovich & Novikov(1971)> for the pressure P and the gravitational mass m, 
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Fig. 3. Gravitational mass Mq versus baryonic mass Mb for some of the stellar configurations shown in Fig|5] The notation is the 
same as in Fig|2](see legend). 



once the EoS P(e) is specified, being e the total energy density and G the gravitational constant. For a chosen central value of the 
energy density e c , the numerical integration of Eqs. Q and @ provides the mass-radius and the mass-central density relations. 
We schematize the entire evolution of a PNS as divided in two main stages. In the first, representing the early stage, the PNS is 
in a hot (T core = 30 -f 40 MeV) stable configuration with a neutrino-trapped core and a high-entropy envelope (s env ^6-^8). The 
second stage represents the end of the short-term cooling where the neutrino-free core possesses a low temperature (T core = 10 
MeV) and the outer part can be considered as a cold crust (BPS+FMT). 

In Fig|5]we report the gravitational mass as a function of the central density (left panel) and stellar radius (right panel) for the 
configurations discussed above. As mentioned before the minimum mass plays a fundamental role in determining the range of 
stability of a PNS (Zeld ovich & Novikov(197T)| > and in particular it is sensitive of the temperature profile assumed. It is found 
that the introduction of an isentropic envelope yields reliable values for the minimum mass for T core = 30 MeV, as shown in 
FigEl For sake of comparison we also show results for a neutrino-free core at the same temperature (dash-double dot line). In 
the neutrino-trapped case the minimum mass rises faster than in the neutrino-free case as the temperature increases. This leads 
to unstable PNS at T core = 40 MeV and s env 8 (dotted line). This result is consistent with the calculations of |Pons et al.( 1999) 
where the temperature profiles do not exceed values of T — 40 MeV. 

The values for the maximum mass lie within a narrow range 1.77 + 1.87Mq for all the considered configurations. A systematic 
collection of results on maximum and minimum mass together with corresponding radii is reported in Tab^ The stellar radii at the 
maximum mass (minimum mass) show an overall trend consisting of an increasing (decreasing) value as T core increases, both for 
the neutrino-free and trapped case. For the second stage the results obtained are very similar to those given in INicotra et al.(2 006) 
for T — 0, 10 MeV and we do not discuss those results further. 

In absence of accretion of matter the baryonic mass is constant during the evolution of the star and we determine a window of 
baryonic mass that can be supported, which ranges from 1.58Mq to 1.98Mq for a PNS with T core = 30 MeV and with trapped 
neutrinos. In Fig[3]we report the gravitational mass Mq versus the baryonic mass Mb in the allowed range. As one can see, for 
a fixed value of Mb the gravitational mass decreases as the PNS evolves to a cold NS. This difference AMq is on average about 
O.IMq and gives a measure of the energy released during the PNS cooling. For pure nucleonic matter metastability 1 does not 
occur. The higher is the temperature of the core the narrower is the window of stability for a PNS, therefore in this stage accretion 
of matter could easily lead to a black hole formation. 

The procedure discussed in section l2T21 works well even for a PNS with hyperons in its interior, giving values for s em , in line with 
those for the pure nucleonic case. Since the development of a finite temperature BHF EoS with interacting hyperons is still in 
progress, we prefer to postpone this discussion to a future work. 




(4) 



1 Here metastability means a decreasing stability as long as the PNS evolves to lower values of temperature. 
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Table 1. Values of maximum and minimum gravitational mass with their corresponding radii and s env for several values of T, 
with and without trapped neutrinos. 



composition 


Tcore 


Senv 


M max /M mm 






[MeV] 


I — ] 




[Km] 









1.859/0.106 


9.7/174.5 


Xy = 


10 




1.867/0.419 


10.0/66.8 




30 


6.13 


1.863/1.320 


12.6/46.3 




40 


7.87 


1.871/1.778 


18.7/47.1 









1.769/0.120 


9.1/61.4 


Y e = 0.4 


10 




1.791/0.825 


10.3/48.9 




30 


6.23 


1.805/1.532 


13.3/39.4 




40 


7.98 


-/- 


-/- 



4. Conclusions 

In this paper we have studied the structure and range of stability of a PNS within a static description. We have taken into account 
two main stages of the evolution characterized by a hot and neutrino-trapped phase and a cold neutrino-free one. The conclusions 
can be summarized as follows. 

1 . A new temperature profile was assumed which resembles those obtained by dynamical calculations for the first few seconds 
of life of a PNS. Particular care has been devoted to the description of the outer part where a temperature drop takes place. 

2. A natural and simple correspondence between temperature of the core and entropy of the envelope has been found. For some 
aspect a hybrid temperature profile between isothermal and isentropic ones has been adopted allowing us to describe the low 
central density part of a PNS, which turns out to be important in determining values of the minimum mass. 

3. The minimum mass rises as the temperature of the core increases. This fact becomes particularly evident when the core is 
neutrino rich, leading to unstable PNS for T core slightly below 40 MeV. 

4. The maximum mass is much more affected by the composition of the core and therefore by the EoS adopted than by the 
global thermal effects. 
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